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ABSTRACT 

We present the near-through mid-infrared flux contribution of thermally-pulsing asymptotic giant branch (TP-AGB) 
and massive red supergiant (RSG) stars to the luminosities of the Large and Small Magellanic Clouds (LMC and 
SMC, respectively). Combined, the peak contribution from these cool evolved stars occurs at ~3-4 /zm, where they 
produce 32% of the SMC light, and 25% of the LMC flux. The TP-AGB star contribution also peaks at ~3-4 /zm 
and amounts to 21% in both galaxies. The contribution from RSG stars peaks at shorter wavelengths, 2.2 /zm, 
where they provide 11% of the SMC flux, and 7% for the LMC. Both TP-AGB and RSG stars are short lived, 
and thus potentially impose a large stochastic scatter on the near-IR derived mass-to-light ( M/L ) ratios of galaxies 
at rest-frame 1-4 /zm. To minimize their impact on stellar mass estimates, one can use the M/L ratio at shorter 
wavelengths (e.g., at 0.8-1 /zm). At longer wavelengths /zm), emission from dust in the interstellar medium 
dominates the flux. In the LMC, which shows strong polycyclic aromatic hydrocarbon (PAH) emission at 8 /zm, 
TP-AGB and RSG contribute less than 4% of the 8 /zm flux. However, 19% of the SMC 8 /zm flux is from evolved 
stars, nearly half of which is produced by the rarest, dustiest, carbon-rich TP-AGB stars. Thus, star formation rates 
of galaxies, based on an 8 /zm flux (e.g., observed-frame 24 /zm at z = 2), may be biased modestly high, especially 
for galaxies with little PAH emission. 
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1. INTRODUCTION 

Stellar masses of galaxies are typically estimated from a 
model mass-to-light (M/L) ratio (e.g., Bell & de Jong 2001), 
and a measurement of the integrated luminosity at optical 
or near-infrared wavelengths. Until recently, galaxy M/L ra- 
tios were thought to be relatively well behaved, especially at 
near-infrared wavelengths where the effects of dust and mas- 
sive main-sequence stars are minimized. However, in the past 
decade, this assumption has been shown to be flawed because of 
the influence of short-lived but extremely luminous thermally- 
pulsing asymptotic giant branch (TP-AGB ) stars (Maraston et al. 
2006; Zibetti et al. 2009, 0.8 < M < 8 M 0 ) and massive red 
supergiant stars (RSG; Dalcanton et al. 2012; Melbourne et al. 
2012). 

Just as the near-IR fluxes of galaxies are often used to estimate 
their stellar masses, mid-IR fluxes can be used to estimate 
their star formation rates. For instance, at rest-frame 8 /zm, 
the flux from starbursts is typically dominated by polycyclic 
aromatic hydrocarbon (PAH) emission which scales with the star 
formation rate (SFR; e.g., Dlaz-Santos et al. 2008). Likewise, 
the rest-frame 24 /zm flux has been shown to scale with the 
IR luminosity of a galaxy and is therefore used as an SFR 
indicator (Chary & Elbaz 2001). However, calibrations of 
these SFR indicators typically ignore the possibility of stellar 
contamination in the flux at these wavelengths. Several works 
have shown that, under certain assumptions, AGB stars could 
contribute significantly to the mid-IR fluxes of galaxies and 
therefore affect the SFRs estimated from their mid-infrared 
fluxes (e.g., Verley et al. 2009; Kelson & Holden 2010; Chisari 
& Kelson 2012). However, as will be discussed in this paper, 
some of their assumptions may be flawed. 


The AGB is a short-lived phase (~1 Myr; e.g., Girardi et al. 
2010) near the end of stellar evolution when a star undergoes 
helium shell burning and swells to a large size, producing a 
relatively cool photosphere and, potentially, a very infrared 
(IR) luminous star. Because of the cool temperatures, complex 
molecules and dust can form in the atmosphere (e.g., Bowen 
1988; Winters et al. 2000, 2003; Schirrmacher et al. 2003; 
Woitke 2006, 2007; Mattsson et al. 2008; van Loon et al. 2008), 
which in turn can drive strong stellar winds (e.g., Sedlmayr & 
Dominik 1995; Elitzur & Ivezic 2001 and references therein) 
that ultimately return as much as half of the stellar mass to 
the interstellar medium. An isolated star will start its AGB 
phase with an oxygen-rich atmosphere (e.g., O-AGB). However, 
during the later stages of AGB evolution, a star undergoes 
a series of thermal-pulses and carbon, dredged up from the 
interior, pollutes the atmosphere (Iben 1983). When the C/O 
ratio exceeds unity, the star is classified as a carbon star (e.g., 
C-AGB), and complex molecules and dust grains form in the 
escaping winds. These dust grains tend to obscure the star at 
optical and near-IR wavelengths, but the same dust emits in 
the mid-IR. Only 20% of the AGB stars in the LMC/SMC are 
estimated to be thermally-pulsing (A. G. Bruzual 2012, private 
communication). However, these TP-AGB stars can be very 
luminous in the infrared, typically exceeding the luminosity of 
the tip of the red giant branch (TRGB), sometimes by several 
magnitudes. 

Red supergiants are helium-burning massive stars that form a 
tight sequence in optical and near-IR color-magnitude diagrams 
(CMDs; e.g.. Figure 1). Stars on this sequence are also termed 
red helium burning stars (RHeB; Dohm-Palmer & Skillman 
2002; McQuinn et al. 2011; Dalcanton et al. 2012). At the 
luminous end, RSGs are truly massive, nine solar masses 
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Figure 1 . J — Ks vs. Ks CMDs for the LMC (gray density plots). Left panel overlays the Padova isochrones (Bertelli et al. 2008; Marigo et al. 2008; Bertelli et al. 
2009; Girardi et al. 2010). Center panel overlays a statistical estimate of the Milky Way foreground stars in the direction of the LMC from TRILEGAL (Girardi et al. 
2005). The right panel shows the Boyer et al. (201 1) classification of RSG (yellow), O-AGB (blue), and C-AGB (red) stars, with selection regions defined by Cioni 
et al. (2006). A luminosity cutoff was made at the TRGB, below which the fluxes of these stars become negligible compared with the RGB. The dust enshrouded 
x-AGB (green) were selected by Boyer et al. (2011) at longer wavelengths; J — Ks >3.1 mag or [3.6] — [8] > 0.8 mag and [8] < 10 mag. Because of their dusty 
envelopes, the x-AGB stars can be fainter than the .ft's -band TRGB and yet be significantly brighter than the TRGB at longer wavelengths. 

(A color version of this figure is available in the online journal.) 



or more, and thus were bom in very recent star formation 
episodes (e.g., <0.1 Gyr). At the luminosity of the K-band 
TRGB, the Padova isochrones predict more modest masses, 
6-7 Mq (Bertelli et al. 2008, 2009; Marigo et al. 2008; Girardi 
et al. 2010), and in fact the RHeB sequence continues to lower 
luminosities for stars of even smaller mass (e.g., McQuinn et al. 
2011). 

Combined, TP-AGB and RSG stars have been shown to 
account for as much 40% of the integrated 1.6 /rm flux, even 
in local galaxies, where the red giant branch (RGB) is well 
developed (Melbourne et al. 2012). In the early universe, when 
the RGB has not had time to form, the TP-AGB and RSG stars 
are expected to contribute even larger fractions (as much as 70%) 
of the near-IR light, even though they represent a tiny fraction of 
the stellar mass (Maraston et al. 2006; Melbourne et al. 2012). 
Thus, these stars must be accounted for to accurately model the 
infrared M/L ratios of galaxies (Zibetti et al. 2009; Ilbert et al. 
2010). 

Not only do these stars contribute to the near-IR luminosities 
of galaxies, but they also can be significant sources of mid-IR 
light (e.g., Boyer et al. 2011). TP-AGB stars, in particular, can 
produce significant amounts of warm dust (e.g., Bowen 1988; 
Schirrmacher et al. 2003). 

These dusty shells both emit in the mid-infrared and obscure 
the TP-AGB stars at optical wavelengths (e.g., Reid 1991). 

The flux contribution of TP-AGB and RSG stars has now 
been well quantified at 0.8 and 1.6 /tm in a set of 23 local 
dwarfs and spirals (Girardi et al. 2010; Melbourne et al. 2012). 
Similarly, the 3.6 and 4.5 /im flux contribution of the most 
extreme (dusty) TP-AGB stars has been quantified in a set of 
six nearby galaxies (Gerke & Kochanek 2013). Likewise, the 
mid-IR (8-24 /urn) flux contribution from carbon stars has been 
estimated from stellar population synthesis models of SINGS 
galaxies (Chisari & Kelson 2012). However, there has not been 
a study that quantifies the flux contributions of these stars across 
the full 1-24 p m wavelength range for a galaxy where complete 
samples of TP-AGB and RSG stars are actually resolved and 
individually measured. 


In this paper, we use complete samples of TP-AGB and 
massive RSG stars in the Magellanic Clouds from the Spitzer 
legacy program “Surveying the Agents of Galaxy Evolution” 
(SAGE; Blum et al. 2006; Meixner et al. 2006; Gordon et al. 
2011; Boyer et al. 2011). We build off of the results of Boyer 
et al. (2011), who showed that the TP-AGB and RSG stars 
contribute significantly to the near- and mid-IR point-source 
fluxes in these galaxies. We now show that their contribution 
to the total (stellar and interstellar) fluxes of the SMC and 
LMC varies significantly from 1 to 24 pm. We also find 
significant differences between our measured flux contributions 
and those from several recent studies in the literature. We find 
a significantly larger TP-AGB contribution at 3.6 pm than is 
found in Gerke & Kochanek (2013) and Meidt et al. (2012a). 
Conversely, we find that the 8-24 pm AGB flux contributions 
estimated in Kelson & Holden (2010) and Chisari & Kelson 
(2012) are likely significantly overestimated. We discuss the 
implications of these results for estimates of stellar mass and 
SFRs of galaxies. 

2. THE DATA 

This paper uses several types of observations: (1) ground- 
based, seeing-limited near-infrared imaging to measure the near- 
infrared fluxes of TP-AGB and RSG stars in the Magellanic 
Clouds, (2) space-based, low-resolution near-infrared imaging 
to determine the total integrated near-infrared luminosities of the 
Magellanic Clouds, and (3) space-based, mid-infrared imaging 
to determine the mid-infrared fluxes of both the stars and the 
integrated galaxy light. In addition to the IR data, we use ground- 
based optical data to explore the flux contribution of short-lived 
massive main-sequence stars at shorter wavelengths. Details of 
these data sets are provided below. All magnitudes are reported 
in the Vega system. 

2.1. Stellar Photometry 

Multi-wavelength stellar photometry of the stars in the Large 
and Small Magellanic Clouds (LMC and SMC, respectively) 
was compiled by the SAGE team. J and K$ (i.e., 1.2 and 
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2.2 fim) near-infrared photometry was obtained from the Two 
Micron All Sky Survey (2MASS; Skrutskie et al. 2006). Mid- 
infrared photometry at 3.6, 4.5, 5.8, and 8 /x m was obtained 
with the Spitzer Space Telescope Infrared Array Camera (IRAC; 
Fazio et al. 2004). Photometry at 24 /im was obtained with 
the Multiband Imaging Photometer for Spitzer (MIPS; Rieke 
et al. 2004). A discussion of the Spitzer observations, image 
reductions, photometry, and source matching for the SMC 
is given in Gordon et al. (2011), and the LMC photometry 
is discussed in Meixner et al. (2006). Optical, BVI, data 
were compiled by the Magellanic Clouds Photometric Survey 
(MCPS; Zaritsky et al. 2002). 

Stellar classifications for stars in the SAGE data were pre- 
sented in Boyer et al. (2011). Here, we briefly summarize the 
method; for a full description of the stellar classification and 
sources of possible contamination, see Boyer et al. (2011). Clas- 
sification was carried out by examining the locations of stars in 
multi-wavelength color-magnitude space. Less-dust-obscured 
TP-AGB and RSG stars were identified in near-infrared CMDs, 
where they are among the most luminous stars, e.g., above the 
AA-band TRGB. At these wavelengths, sequences of carbon- 
rich and oxygen-rich AGB stars (C-AGB and O-AGB stars, 
respectively) separate by J — K$ color, with C-AGB stars typi- 
cally being redder than O-AGB stars. Likewise, RSG stars form 
a tight sequence slightly blueward of the O-AGB stars. These 
sequences are clearly visible in the J — Ks CMD of the LMC 
shown in Figure 1 which also shows the color-magnitude bound- 
aries used to identify these different sequences (adopted from 
Cioni et al. 2006). While these boundaries are useful for identi- 
fying different classes of stars, we note that there will be some 
overlap between the populations. Boyer et al. (201 1) provide a 
detailed discussion of classifications including estimates of the 
numbers of misclassified stars, which are small. We note that 
the anomalous O-rich AGB stars reported by Boyer et al. (2011) 
are included with the O-AGB sample here. 

The most dust enshrouded TP-AGB stars are significantly 
extinguished even in the near-IR, and therefore they must be 
identified at longer wavelengths (e.g., Reid 1991). Boyer et al. 
(2009) found that short wavelength searches can miss up to 
40% of the TP-AGB stars. These so-called extreme AGB stars 
(x-AGB) were identified by their colors at mid-IR wavelengths 
(3.6-8 /rm). Most x-AGB stars are C-rich (van Loon et al. 
2005). Combining these two identification methods allows for 
the selection of near-complete samples of TP-AGB stars (e.g., 
these classifications do not miss the most dust enshrouded stars, 
which can be missed by classifications that rely entirely on 
shorter wavelengths). 

While these techniques will recover the bulk of the luminous 
TP-AGB stars, some TP-AGB stars will be missed because 
they are at a minimum in their thermal pulsation cycle and 
are fainter than the TRGB. Models of the TP-AGB stars in 
the Magellanic Clouds based on the Padova stellar evolution 
tracks reveal that <10% of them lie below the TRGB due to 
the thermal pulse cycle (e.g., Marigo & Girardi 2007). We have 
included these subluminous TP-AGB stars in the uncertainties 
of the integrated TP-AGB fluxes by assuming that they have 
the luminosity distribution of the top magnitude of the RGB. 
We added this uncertainty in quadrature to the other sources 
of uncertainty, increasing the upper error bars in Table 1 and 
Figure 2 by <3%. 

We also note that there are AGB stars that are not thermally- 
pulsing. These are lower-luminosity AGB stars and also lie 
below the TRGB. These stars cannot be easily selected in CMD 


space because they overlap the RGB stars. However, as these 
AGB stars are also relatively rare compared with the equally 
luminous RGB, they will not contribute significantly to the total 
flux of the galaxy. 

The total numbers of SMC and LMC stars within each stellar 
classification bin are presented in Table 1 . 

2.2. Foreground Contamination 

The CMDs of the LMC and SMC also contain some Milky 
Way foreground stars with the colors and luminosities of 
Magellanic Cloud AGB and RSG stars. We applied a statistical 
correction for foreground contamination of Milky Way stars at 
each wavelength. We used TRILEGAL (Girardi et al. 2005) 
to model the numbers and fluxes of Milky Way stars in the 
direction of the Magellanic Clouds in the same filters as our 
data. The TRILEGAL photometry was then run through the 
same classification scheme as the RSG and TP-AGB stars. The 
contamination is largest for the RSG population of the SMC 
(~25%). In the LMC, the RSG contamination is negligible 
because the higher metallicity of the LMC causes all of the RSG 
and AGB stars to be redder than they are in the SMC, resulting 
in less overlap with the foreground sources in the near-IR CMD 
(see Figure 1). Likewise, the AGB populations in both galaxies 
are effectively foreground contamination free. The estimated 
numbers of foreground stars within each stellar class are also 
presented in Table 1 . 

Note: stars in the Milky Way foreground clusters 47 Tuc and 
NGC 362, in the direction of the SMC, were removed from the 
data (see Boyer et al. 201 1). 

2.3. Integrated Galaxy Photometry 

Measuring the total luminosity of the Magellanic Clouds is 
challenging because of their very large areas on the sky. Ground- 
based, near-IR images are dominated by thermal background 
flux, and thus accurately discerning sky from unresolved galaxy 
becomes nearly impossible. Fortunately, these issues are largely 
eliminated by going to space, where the IR thermal background 
is minimal and COBE/ DIRBE (Hauser et al. 1998) provides a 
large beam on the sky (nearly a degree). 

We use the COBE I DIRBE integrated flux density estimates 
of Israel et al. (2010) for the Magellanic clouds at 1.25, 2.2, 
and 3.5 /rm. The SMC flux is measured within a 12?3 x 12?3 
area, and the LMC flux covers a 15° 3 x 15? 3 area. These 
measurements were carefully corrected for foreground (and 
background) contamination by Israel et al. (2010). The Israel 
et al. flux density estimates for the Clouds are presented here in 
Table 1. 

Similarly, the global SMC and LMC IR fluxes from the 
Spitzer data at 3.6, 4.5, 5.8, 8, and 24 /rm were made by Gordon 
et al. (2011) and K. D. Gordon et al. (in preparation). These 
fluxes were measured within circular apertures of radius 2? 25 
for the SMC and 3.75 for the LMC. For both galaxies, a sky 
annulus with a 0?2 width was used to estimate and subtract 
the background/foreground. To within the uncertainties, the 
integrated DIRBE fluxes at 3.5 /xm from Israel et al. match the 
Spitzer fluxes at 3.6 /xm from Gordon et al. The Spitzer fluxes 
are provided in Table 1 . 

The spectral energy distributions (SEDs) of the SMC and 
LMC are shown in Figure 2. Also overlaid are dust model 
predictions from Bot et al. (2010). While the stellar flux is 
declining with wavelength, the warm dust and PAH components 
are increasing. 
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Table 1 

Integrated Flux Density of the SMC and LMC from 1 to 24 (tm and Percentage of that Flux Contributed by TP-AGB and RSG Stars 


Type 

Stars 3 

Fgnd b 

1 .2 pm c 

2.2 pm c 

3.5 /im d 

3.6 fi m e 

4.5 /xm e 

5.8 /xm e 

8 /xm e 

24 p.m c 

SMC 

















Integrated Flux (Jy) 







670 ± 80 f 

525 ± 65 f 

280 ± 40 f 

300 ± 208 

220 ± 108 

220 ± 10 s 

200 ± 10 g 

350 ± 108 







Fractional Contribution (%) 



O-AGB 

3 353 

4 

r 7 + 1.1 

J *'-0.7 

O 4+1.5 

8.02‘ifi 

7 5 +l1 
'• :, -o.5 

6 O +0 - 8 
o>u -0.3 

4 5 + 0.6 
— 0.2 

O 1+0.4 

Q 14 +O.OI 
u - A ^-o.oi 

C-AGB 

1 597 

<1 

3-02°o 4 4 

6 6 +0 - 9 
0.8 

9.12V 4 3 

Q C+0.8 

“•^-o.e 

7 i + 0.6 

' ' 1 —0.3 

5 3 +0 - 4 
J - 0.2 

4 8 +0-4 
^•°- 0.2 

0 31 +0 02 
u - :)A -o.oi 

x-AGB 

309 

<1 

°- 2 -oi 

1 2 +01 
A-Z - 0.1 

4 8 +0 - 7 

4 5 + 0.3 

6 5 +0 - 3 
°- J -0.3 

6 9 +0 - 3 
°- y -0.3 

7 6 +0 - 4 
/ -°-0.4 

1 50+ 004 
1 - r,u -0.04 

RSG 

2 340 

598 

8.97V 1 ! 

10-97V 4 4 

10 . 02 V 4 

9.42 ° 0 6 6 

n ft+ 0.3 
' ■ — 0.3 

5 < 7 + 0.3 

o 7 + 0.2 

3 -'-0.2 

Q 00 +0.01 

U~3.3_0.01 

Total 

7 599 

602 

n.s/VI 

27. L 2 ! ,9 

31.92V 4 3 

29.82V 5 o 

27 2 +2-3 
z/ - z - 0.6 

77 o+ 2. 0 

zz> — 0.5 

19 3 +21 
— 0.5 

2 l + 0- 2 
Z - A — 0.1 

LMC 

















Integrated Flux (Jy) 







4520 ± 650 f 

3 770 ± 540 f 

2 190 ± 300’ 

f 2 080± 100 h 

1 350 ± 70 h 

1 850 ± 200 h 

5 980 ± 300 h 

7 640 ± 200 h 







Fractional Contribution (%) 



O-AGB 15 635 

<1 

6.87 2 i 6 o 

9.92V 1 ! 

S-S-d.s 

Q 7+3.2 
^-O.S 

O 0+3.1 

o.D -0 4 

4 6 +l7 
^■°-0.5 

0 9+0-3 
u - y -0.1 

q 1 1 +0.22 
U - A A -0.01 

C-AGB 

5 628 

<1 

2 6 +0 - 4 
z -°-0.4 

6 0 +0 - 7 
— 0.6 

2 -67°o 7 7 

o i +0.5 

o^-O^ 

7 O+0.4 
/ - :> -0.4 

o 7+0.4 
—0.4 

1 Q+0. 1 

A - u -o.i 

0 1 1 + 0- 01 
u - AA -o.oi 

x-AGB 

989 

<1 

0 l +01 
u - A -o.i 

0 7 +01 
u>/ -0.1 

2 Q +0 - 3 
Z - y -0.3 

3 0 +0 - 2 
•^-o.i 

5 2 +0 ’ 3 
J ' z -0.3 

4 5+0.5 

^-o.s 

1 4 + 0.1 
A *^- 0.1 

0 55 + 0.01 

u.DJ-ooi 

RSG 

4 361 

13 

5 5 +0 - 8 
J - 0.8 

6 9 +0 ’ 7 
— 0.7 

6 0 +0 - 5 
°- u — 0.5 

6-32°o 3 3 

5 5 +0 - 3 
— 0.3 

3 0 +0 - 3 
•^-o.s 

0 6 +01 
u -°-o.i 

n 1 A+0.01 

Total 

26 613 

13 

! 5.12V 3 

23.52V 4 

25.22V V 

26.62 3 o 3 7 

26 3 +3 - 2 
zo -0.7 

1 c zr+1.9 
15.6_o9 

o q+ 0.4 
■^-o.i 

0 7 +0 - 2 


Notes. 

a Number of stars of given type detected at 3.6 /cm; the number varies slightly with wavelength. 
b Number of Milky Way foreground stars as estimated by TRILEGAL (Girardi et al. 2005). 
c Integrated Flux from COBE/DIRBE (at 1.25 and 2.2 (im), Stellar Flux from 2MASS (at 1.2 and 2.2 //m). 
d Integrated Flux from COBE/DIRBE (at 3.5 // in). Stellar Flux from Spitzer (at 3.6 gm). 
e Integrated Flux from Spitzer , Stellar Flux from Spitzer. 
f From Israel et al. (2010). 

8 From Gordon et al. (201 1 ). 
h From K. D. Gordon et al. (in preparation). 


We compare these integrated galaxy fluxes gathered from 
the literature to the fluxes from the TP-AGB and RSG stars as 
measured from 2MASS and Spitzer photometry. The 2MASS J 
and K s bands are similar to the COBE / DIRBE 1.25 and 2.2 pm 
filters, however, the DIRBE filters are 0.05 pm (15%) wider at 
the long wavelength end. Similarly, the CDBZs/DIRBE 3.5 pm 
filter is slightly bluer and about 25% wider than the Spitzer 
3.6 pm filter. These slight filter mismatches may increase the 
uncertainties on the measured AGB and RSG flux fractions 
by modest amounts. There are no filter mismatches for the 
3.6-24 pm flux ratios, as all measurements are made from the 
same Spitzer images. 

2.4. Uncertainties 

In the following section, we report the fractional flux contri- 
bution of TP-AGB and RHeB stars to the Magellanic Clouds as 
given by: 

Stellar Flux — Foreground Stellar Flux 

Flux Fraction = 5 . (1) 

Integrated Galaxy Flux 

The uncertainties for these measurements are constructed 
from the following components. First, we assume the published 
uncertainties for the integrated fluxes of the Magellanic Clouds 
at each wavelength. Second, we compile the photometric un- 
certainty in the integrated stellar flux from the photometric 
uncertainty of each star added in quadrature. Third, we esti- 
mate the uncertainty on the foreground correction by taking 
the Poisson uncertainty on the numbers of foreground stars 


in each classification multiplied by the typical flux of those 
stars. We propagate these uncertainties through Equation 1. 
Fourth, we include an estimate of the uncertainty introduced 
by AGB stars that are fainter than the TRGB as described in 
Section 2.1. 

Finally, we have quantified how much of the TP-AGB and 
RSG population we have missed due to limited spatial coverage 
within the apertures we used to measure each galaxy’s inte- 
grated flux. In both galaxies, the TP-AGB and RSG stars follow 
a spherical distribution centered on the Bars (e.g., Blum et al. 
2006; Boyer et al. 2011). For the SMC, the SAGE footprint in- 
cludes an additional strip of coverage that extends >6 deg along 
the SMC Tail, and for the LMC, the SAGE coverage stretches 
along the Disk, >4.5 deg from the Bar’s center. We can use this 
extended coverage to measure the radial profile of AGB stars 
in the Bars by counting the number of AGB stars in concentric 
ellipses centered on the Bars. We find that the AGB population 
is rather tightly concentrated, dropping to very small numbers 
by about 2 deg (along the minor axes) from the population 
centers. By extrapolating these radial profiles in all directions, 
we estimate that we have missed about 500 SMC AGB stars 
and 900 LMC AGB stars (10 SMC RSGs and 20 LMC RSGs) 
by restricting our apertures to 2.25 and 3.75 deg for the SMC 
and LMC, respectively. Assuming that the missed stars have a 
flux distribution that is similar to the measured TP-AGB /RSG 
population, this amounts to an extra <1% uncertainty on the 
integrated AGB (and RSG) flux, which we have included in the 
upper uncertainties in Table 1 and Figure 2. 
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Figure 2. Fractional contribution from TP-AGB stars (open symbols) and RSG stars (plus symbols) to the integrated near- and mid-IR fluxes of the SMC and LMC 
(top and middle panels). The contributions divided by subspecies are given in the middle panels, while the total contributions from all TP-AGB stars (circles) and 
TP-AGB + RSG (asterisk) are shown in the top panels. The bottom panel shows the integrated SEDs of the LMC and SMC and a model of the dust emission from 
Bot et al. (2010). While the stellar flux declines with wavelength (dotted line), the dust emission increases. The flux contribution of TP-AGB stars peaks at ~3-4 /im. 
The contribution from RSG stars peaks at shorter wavelengths (roughly 2.2 /xm). Combined, these evolved stars contribute 32% of the SMC light at 3-4 /xm and 19% 
at 8 /xm. The flux contributions of the TP-AGB and RSG at 24 /xm are small, <3%, in both galaxies. 

(A color version of this figure is available in the online journal.) 


3. RESULTS: THE MULTI- WAVELENGTH 
CONTRIBUTION OF TP-AGB AND RSG STARS TO THE 

INFRARED FLUX OF THE MAGELLANIC CLOUDS 

Figure 2 shows the fraction of the integrated galaxy light 
produced by TP-AGB (open symbols) and RSG (plus symbols) 
stars in the SMC and LMC as a function of wavelength from 1 to 
24 /xm. These results are also presented in Table 1. At 3.6 /xm, 
these rare stars contribute over 30% of the SMC flux, even 
though they represent a negligible fraction of the stellar mass 
of the galaxy (numbers of stars of each type are also shown in 
Table 1). In the LMC, they contribute more than a quarter of 
the 3.6 /xm light. At 24 /xm, the stellar contribution to the total 
luminosity is <3% for both galaxies. 

The flux contribution from RSG stars, which are bluer than 
AGB stars, peaks near ~2 /xm. In the SMC, they are responsible 
for over 10% of the L-band light. In the LMC, the fraction 
is lower, 7%. Part of this difference is no doubt due to the 
detailed star formation histories of these two galaxies over 
the last 0.1 Gyr, with the SMC having a higher specific star 
formation rate over that period (e.g., Harris & Zaritsky 2004, 
2009). 

The flux contribution from TP-AGB stars peaks in the Spitzer 
IRAC bands between 3 and 4 /xm, reaching 15%-20% for both 
the LMC and SMC. The similarity of the results for the two 


galaxies suggests that over the longer timescale of 0. 1-2 Gyr, 
these galaxies have had similar star formation histories, as has 
been inferred by studies of their CMDs (Harris & Zaritsky 2004, 
2009). 

Figure 2 also shows the TP-AGB populations subdivided into 
the three categories previously discussed: O-AGB, C-AGB, and 
x-AGB. Despite the fact that the O-AGB population outnumbers 
the C-AGB population by a factor of 1.8 in the SMC and 
2.4 in the LMC, the contribution to the global flux of the 
O-AGB and C-AGB populations are similar at X > 3 /xm, with 
each contributing roughly equal parts to the total TP-AGB flux. 
However, the contribution from O-AGB stars peaks at shorter 
wavelengths (near 2.2 /xm) than the C-AGB star contribution 
(near 3.6 /xm). In contrast, the x-AGB stars account for less 
than 1% of the global flux at X < 2/xm, but exceed the flux 
contribution of the other two TP-AGB types at X > 4 /xm. This 
is true even though the x-AGB population is numerically tiny, 
comprising <4% of the total number of TP-AGB stars in both 
galaxies. 

4. DISCUSSION 

We present, for the first time, the near- through mid-IR flux 
contributions of complete samples of TP-AGB and RSG stars 
to the integrated luminosities of the Magellanic Clouds. We find 
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that these stars can contribute non-negligible flux fractions from 
1 to8 jim (e.g., 15%-30%),and smaller fractions at 24 ;xm (e.g., 
<3%). Below, we compare these results to previous work and 
discuss the implications for estimating stellar masses and star 
formation rates of galaxies. 

4.1. The Importance ofTP-AGB and RSG Stars 
at Near-IR Wavelengths 

Rest-frame near-IR fluxes typically anchor the stellar mass 
estimates of galaxies (e.g., Maraston et al. 2006; Conroy 
et al. 2009; Zibetti et al. 2009; Eskew et al. 2012). At these 
wavelengths, dust obscuration is significantly reduced compared 
with optical light, and a galaxy’s near-IR luminosity is much 
less affected by the short-lived O- and B-type stars produced by 
current star formation. Instead, the near-IR is typically thought 
to be dominated by older, easily modeled RGB stars. 

However, we have shown that TP-AGB and massive RSG 
stars can also account for large fractions of the near-IR light 
(e.g., >30% at 3.6 /.mi for the SMC), even though they contribute 
essentially nothing to the total stellar masses of galaxies. A 
galaxy with the same 3.6 jim luminosity as the SMC but 
comprised of entirely old stellar populations (e.g., >3 Gyr) must 
be nearly twice as massive as the SMC to make up for the 
missing stellar luminosity provided by the TP-AGB and RSG at 
that wavelength. Thus, a galaxy’s flux at 3-4 jim (and hence its 
M/L ratio) is highly sensitive to the star formation history. 

The large TP-AGB and RSG near-IR flux contributions we 
measure for the Magellanic Clouds fit well with the contribu- 
tions reported in Melbourne et al. (2012) for a sample of 23 
nearby dwarfs and small spirals, which ranged from 5% to 40% 
at 1.6 jim. Similarly, our results are well matched to the pre- 
dictions of Bruzual (2011), which discussed the contribution of 
AGB stars to the K-band flux of their stellar population synthesis 
models. Their results for a galaxy with Z = 0.008 and a de- 
clining SFR with an e-folding time of 1 Gyr are consistent with 
what we see in the Magellanic Clouds. These large flux contri- 
butions suggest that galaxies with even modest populations of 
TP-AGB and RSG stars will have very different near-IR M/L 
ratios than galaxies lacking these populations, making near-IR 
based stellar mass estimates difficult. 

However, some authors find significantly smaller near-IR 
flux contributions from TP-AGB stars. For instance, Gerke & 
Kochanek (2013) calculate a smaller 3.6 pm flux contribution 
from the TP-AGB in the Magellanic clouds, only 5% of the 
integrated light. They selected AGB stars using the [3.6] — [4.5] 
color, however, they selected a CMD region that is comprised 
predominantly of x-AGB stars. Thus, they miss the bulk of 
the TP-AGB altogether. The top panel of Figure 3 shows the 
distribution of [3.6] — [4.5] colors of TP-AGB stars in the 
LMC. The Gerke & Kochanek (2013) sample has a mean color 
similar to the x-AGB alone (green dot-dashed line), which are 
significantly redder than the overall sample (dotted line). If we 
compare only the x-AGB stars in our sample to their result, then 
we find a similarly small contribution of ~3%. Including the 
full TP-AGB sample results in a much higher flux contribution 
at 3.6 jim, as discussed above. 

Meidt et al. (2012a) show a similar result to Gerke & 
Kochanek (2013), only they used unresolved Spitzer 3.6 and 
4.5 pm observations of galaxies in the Spitzer Survey of Stellar 
Structure in Galaxies (S 4 G; Sheth et al. 2010). Their work 
spatially separates (on a pixel-by-pixel basis) older stellar 
populations from warm dust, PAH emission, and evolved stars, 
based on the colors of regions in the Spitzer data. Thus, they 
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Figure 3. Histograms of the [3.6] — [4.5] (top) and Ks — [8] (bottom) colors of 
the AGB (black) and RSG (yellow) stars in the LMC. Also shown are various 
color selections used in the literature to select AGB populations. The top panel 
shows the [3.6] — [4.5] color range used by Meidt et al. (2012b; shaded region) 
to isolate AGB and RSG stars, which roughly matches the mean selection used 
in Gerke & Kochanek (2013). Both of these selections are significantly redder 
than the mean AGB and RSG colors in the LMC (dotted line; and the SMC, see 
Boyer et al. 2011). The lower panel shows the typical carbon star color used in 
the Kelson & Holden (2010) models (dashed line), which is significantly redder 
than the mean color for the x-AGB stars (green dash-dotted line) and over 4 
mag redder than the mean carbon star color (red triple-dot-dashed line). 

(A color version of this figure is available in the online journal.) 

are able to nicely constrain the warm dust and PAH emission to 
roughly 5%— 15% of the 3.6 pm flux in a sample of six nearby 
galaxies. Likewise, they find that pixels dominated by evolved 
AGB and RSG stars account for only another roughly 5% of the 
3.6 jim light. However, as with the Gerke & Kochanek (2013) 
result, their assumed AGB color of [3.6] — [4.5] = 0.24-0.76 
(gray shaded region of Figure 3) is significantly redder than the 
bulk of the TP-AGB stars in the LMC which have a mean color 
of [3.6] — [4.5] = 0.0 (red triple-dot-dashed line). In fact, the 
mean [3.6] — [4.5] color for TP-AGB and RSG in the LMC is 
[3.6] — [4.5] = —0.07 (dotted line), exactly the color Meidt et al. 
(2012a) assign to old stars. Unfortunately, that means that Meidt 
et al. (2012a) have only selected the reddest TP-AGB stars and, 
as a result, they are likely underestimating the contribution of 
TP-AGB stars at 3.6 jim. 

The results we have presented here suggest that the near-IR 
may not be optimized for stellar mass estimates of galaxies. 
The 3.6 jim band, in particular, is the most problematic for the 
effects of short-lived TP-AGB and RSG stars, which can impose 
a stochastic variability on the near-IR M/L ratio of stellar 
populations (e.g., Melbourne et al. 2012). However, it should be 
possible to minimize the effects of the TP-AGB and RSG stars 
on estimates of the near-IR M/L ratios of galaxies by using bluer 
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Figure 4. Contribution of luminous stars (with fluxes > than the TRGB) to the 
total point-source flux (measured to 2.5 mag below the TRGB) as a function 
of wavelength. At the blue end the O and B stars contribute significantly to the 
point-source flux. At the red end, RSG and TP-AGB stars contribute significantly 
to the total. As all of these stars are short lived, they will result in a rapidly 
changing M/L ratio at these wavelengths. The impact of luminous short-lived 
stars on the stellar mass estimates of galaxies can be minimized by using the 
M/L ratio at 0. 8-1.0 fim, where this plot shows a minimum. 



wavebands. The TP-AGB and RSGs flux contribution at 1 pm 
is roughly half that of the 3.6-/xm flux contribution. Of course, 
bluer wavelengths will be more impacted by extinction from 
dust and emission from hot O and B stars, but the I or J bands 
may prove to be the ideal compromise. To test this hypothesis, 
we compare the point- source flux above the TRGB to the total 
point-source flux (which is measurable to 2.5 mag below the 
TRGB) as a function of wavelength from 0.5 to 3 .6 pm. Figure 4 
displays the results, and shows that there is indeed a minimum 
flux contribution from luminous short-lived stars as a function 
of wavelength. The minimum occurs between 0.8 and 1 /xm. 
At shorter wavelengths, massive main-sequence O and B stars 
contribute large fractions of light. At longer wavelengths, the 
TP-AGB and RSG stars contribute large fractions. This plot 
suggests that rest-frame / through J bands will indeed provide 
M/L ratios that are the most stable against stellar population 
history and the influence of these short-lived luminous stars. 

4.2. The Importance of TP-AGB and RSG Stars 
at Mid-IR Wavelengths 

Because mid-IR light is produced by warm dust and PAH 
emission in star forming regions, mid-IR luminosities of galax- 
ies are often used as proxies for their current star formation rates. 
However, we have shown that the TP-AGB and RSG stars can 
also contribute to the mid-IR fluxes of galaxies, with the bulk 
of the mid-IR contribution coming from the rare but extremely 
dusty x-AGB stars. In the SMC, the contribution of TP-AGB 
and RSG stars at 8 /xm is »H9%! If these x-AGB stars are ne- 
glected in the accounting of the mid-IR light, then an SMC SFR 
based on 8 /xm light alone will be biased high. 

This result has impacts on the interpretation of thousands of 
mid-IR based star formation rates, for example at z = 2, Spitzer 


24 /xm imaging traces rest-frame 8 /xm, where the TP-AGB 
is potentially an important flux contributor. However, the size 
of the effect may be small in practice. For instance, in the LMC, 
the TP-AGB + RSG 8 /xm contribution is <5%. Compared with 
the SMC, the LMC has a significantly larger contribution from 
warm dust and 8 /xm PAH emission (Bot et al. 2010, and Figure 2 
here). Samples of the extremely dusty z — 2 galaxies detected 
by Spitzer have been shown to have large PAH flux contributions 
(e.g., Desai et al. 2009), so the correction for TP-AGB stars may 
be small in the high-- 7 luminous infrared galaxies. However, 
it is possible that TP-AGB stars artificially elevate the rest- 
frame 8 /xm estimated SFRs of galaxies, especially for lower- 
metallicity systems like the SMC, where (1) PAH emission 
is limited, (2) C-AGB stars are more common, and (3) star 
formation rates are likely smaller. 

We are not the first to point out that the TP-AGB may 
be important at mid-IR wavelengths. For instance, Kelson & 
Holden (2010) and Chisari & Kelson (2012) recently proposed 
TP-AGB stars as the primary source of mid-IR flux in galaxies, 
and may account for the rise of luminous infrared galaxies 
with redshift. Similarly, Verley et al. (2009) suggested that 
dusty shells surrounding AGB stars could account for the bulk 
of diffuse 24 /xm emission in the outskirts of M33. While 
our results show the importance of TP-AGB stars at mid-IR 
wavelengths, these earlier works likely overemphasize the TP- 
AGB contribution, especially at 24 /xm, where we show that 
TP-AGB stars contribute less than 3% of the Magellanic Cloud 
flux. 

The large mid-IR TP-AGB contributions predicted by these 
previous works (e.g., >50% of the 8 and 24 pm flux) likely 
arose from their adopted mid-IR colors and fluxes of C-AGB 
stars. Kelson & Holden (2010) and Chisari & Kelson (2012) 
assume that C-rich AGB stars have an average K — [8.8] 
color of 5.66 ±1.16 mag, derived from Galactic carbon stars 
observed in the mid-IR by Guandalini et al. (2006). Verley et al. 
(2009) assume an average 24 pm flux derived from carbon stars 
observed with Spitzer in the Magellanic Clouds by Groenewegen 
et al. (2007). However, both the Guandalini et al. (2006) and 
Groenewegen et al. (2007) samples are heavily biased toward 
the reddest x-AGB stars and completely exclude the bulk of 
the C-AGB population, which is much bluer. This issue is 
demonstrated in the bottom panel of Figure 3, which plots 
histograms of the K$ — [8] color of the TP-AGB stars in the 
LMC. The color adopted by the Kelson & Holden (2010) models 
is shown as a dashed line and is significantly redder than the 
bulk of the LMC C-AGB stars. When including the complete 
samples of C-rich AGB stars (C-AGB plus x-AGB, red and 
green distributions, respectively) from the SAGE survey, we 
find ( Ks — [8]) x+ c = 1.38 mag in the SMC and 1.31 mag 
in the LMC (triple-dot-dashed line in Figure 3). This is over 
4 mag bluer than the Kelson & Holden models (also used 
in Chisari & Kelson 2012). A similarly sized discrepancy is 
obtained when comparing the SAGE AGB sample to the adopted 
24 pm flux for AGB stars in Verley et al. These revised colors 
would result in a much smaller contribution of TP-AGB stars at 
mid-IR wavelengths, as is seen in Figure 2 for the Magellanic 
Clouds. 

While some x-AGB stars in the SAGE Magellanic Cloud 
sample do show colors as red as those used in the Kelson 
& Holden models, together they have an average color that 
is still 2 mag bluer than the models, (Ks — [8]} x , smc = 
3.25 mag and (Ks — [8]) x , lmc = 3.41 mag (dot-dashed line 
in top panel of Figure 3). In addition, these x-AGB stars 
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contribute almost nothing to the K-band light of the Magellanic 
clouds (1% or less). The bluer O-AGB and C-AGB stars 
produce roughly 20 times the K - band light of the x-AGB stars. 
The Kelson & Holden models rely on the K-band flux to 
estimate the AGB contribution at longer wavelengths, and they 
assume that the x-AGB stars are contributing nearly 50% of the 
K-band light. Because of these assumptions, they are likely 
significantly overestimating the AGB flux contribution at longer 
wavelengths. 

5. CONCLUSIONS 

The contribution of TP- AGB and RSG stars to the integrated 
fluxes of the Magellanic Clouds varies with wavelength. We 
find that their combined contribution peaks between 3-4 /x m, 
where they contribute 32% of the SMC flux and ~25% of 
the LMC flux. At 1.2 /xm, the contributions are smaller, ~ 1 8% 
and ~ 1 5% for the SMC and LMC. At longer wavelengths, the 
contributions also decline, but still constitute 19% of the 8 /xm 
light of the SMC (only 4% for the LMC where warm dust and 
PAH emission are stronger). The bulk of this 8 jim stellar flux is 
from the most dust-enshrouded AGB stars, the so-called x-AGB 
stars, which are predominantly carbon stars near the end of the 
AGB phase. At 24 /xm, the AGB star contribution to the total 
flux is less than 3% for each galaxy. 

The flux contributions of the TP-AGB and RSG are significant 
at 3.6 /x m. Thus, these stars have an impact on the M/L ratios of 
the Magellanic Clouds at near-IR wavelengths. Because these 
stars are short-lived, we can expect that the near-IR M/L ratios 
will also vary significantly over time. To minimize the impact 
of short lived TP-AGB and RSG stars on stellar mass estimates 
of galaxies, shorter wavelength observations are preferred. 
However, at optical wavelengths, short-lived, superluminous O 
and B stars play a similar role as the AGB and RSG stars in the 
near-IR. We find that the optimal wavelengths for minimizing 
the impact of luminous short-lived stars on the M/L ratio of 
galaxies is between 0.8 and 1 /xm. 

The non-negligible flux from the TP-AGB at 8 /x m suggests 
that star formation rates based on rest-frame 8 /xm flux (e.g., 
Spitzer 24 /xm imaging at z = 2) may be biased slightly high. 
However, this effect is most important for systems with little 
PAH emission (like the SMC). For systems with significant 
warm dust and PAH emission, the TP-AGB contribution at 8 /xm 
appears to be small, and thus for the LMC the contribution 
is <4%. 

Special thanks to Karl Gordon for providing integrated fluxes 
of the LMC in the Spitzer bands, and to Caroline Bot for shar- 
ing the SED and model fits of the Magellanic Clouds. Thanks to 
Gustavo Bruzual for sharing updated single stellar population 
models of the Magellanic Clouds. We also thank Barry Madore, 
Mark Seibert, and Wendy Freedman for inspiring efforts to study 
the impact of AGB stars on mid-infrared galaxy fluxes. This 
work is based in part on observations made with the Spitzer 
Space Telescope , which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology under a contract with 
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